Centrality dependence of charged-neutral particle fluctuations in 158 A GeV 

208pb+208pb collisions 
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Results on the study of localized fluctuations in the multiplicity of charged particles and photons 
produced in 158-y4 GeV/c Pb+Pb collisions are presented for varying centrality. The charged 
versus neutral particle multiplicity correlations in common phase space regions of varying azimuthal 
size are analyzed by two different methods. Various types of mixed events are constructed to 
probe fluctuations arising from different sources. The measured results are compared to those 
from simulations and from mixed events. The comparison indicates the presence of non-statistical 
fluctuations in both the charged particle and photon multiplicities in limited azimuthal regions. 
However, no correlated charged-neutral fluctuations, a possible signature of formation of disoriented 
chiral condensates, are observed. An upper limit on the production of disoriented chiral condensates 
is set. 

PACS numbers: 25.75.+r,13.40.-f,24.90.-|-p 

I. INTRODUCTION Synchrotron (SPS) and the Relativistic Heavy-Ion 



The large number of particles produced in rel- 
ativistic heavy-ion collisions at the Super Proton 
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Collider (RHIC) provide an opportunity to analyze 
and study, on an event-by-event basis, fluctuations 
in physical observables, such as particle multiplic- 
ities, transverse momenta, and their correlations. 
These studies provide information on the dynam- 
ics of multi-particle production and may help to re- 
veal the phase transition from hadronic matter to 
quark-gluon plasma (QGP) P, The formation 
of hot and dense matter in high energy heavy-ion 
collisions also has the possibility to create matter in 
a chiral symmetry restored phase in the laboratory. 
After the initial stage of the collision, the system 
cools and expands returning to the normal QCD vac- 
uum in which chiral symmetry is spontaneously bro- 
ken. During this process, a metastable state may be 
formed in which the chiral condensate is disoriented 
from the true vacuum direction. This transient state 
would subsequently decay by emitting pions coher- 
ently within finite sub volumes or domains of the 
collision region. The possibility of formation of dis- 
oriented chiral condensate (DCC)has been discussed 
extensively in recent years IIIIItIIIII llollTH . 
The detection and study of the DCC state would 
provide valuable information about the chiral phase 
transition and the vacuum structure of strong inter- 
actions. 

Theoretical studies 0, IE suggest that isospin 
fluctuations caused by formation of a DCC would 
produce clusters of coherent pions in localized phase 
space domains. The formation of DCC domains 
would be associated with large event-by-event fluc- 
tuations in the ratio of neutral to charged pions. The 
probability distribution of the neutral pion fraction, 
/, in such a DCC domain has been shown ^ to fol- 
low the relation: 

P{f)^^ where f^N^o/N^ (1) 

which is quite different from that of the normal pion 
production mechanism. For the normal pion produc- 
tion, where the production of 7r'^,7r+, tt^ are equally 
probable, the / distribution is binomial peaking at 
1/3. 

In the experimental search for localized domains of 
DCC a practical approach is to search for events with 
large and localized fluctuations (localized in pseudo- 
rapidity (77) and azimuthal angle (0) ) in the ratio 
of the number of photons to charged particles which 
would directly reflect fluctuations in the neutral to 
charged pion ratio. Typical event structures would 
be similar to the Centauro and anti-Centauro events 
reported by the JACEE collaboration Results 
from other cosmic ray experiments have not ruled 
out the possibility of the DCC formation mechanism 
p^sf. The accelerator based studies carried out in 
p — p IT3| and heavy ion '15|, Il6| reactions have in- 
vestigated particle production over extended regions 



of phase space. These analysis were not sensitive to 
the presence of small domains of DCC localized in 
phase space. A first search for evidence of localized 
domains of DCC has been carried out at the SPS 
by the WA98 experiment in a detailed study of cen- 
tral Pb-f-Pb events The analysis showed the 
presence of localized non-statistical fluctuations in 
the multiplicity of both photons and charged parti- 
cles. However, the charged-neutral fluctuations were 
found not to be correlated event- by-event, as would 
be expected for a DCC production mechanism. An 
upper limit on the frequency of DCC formation in 
central Pb-|-Pb collisions was set. Recently there 
have been theoretical suggestions to look for DCC 
formation in events for intermediate centralities [I^l ■ 
In this paper we present first results on the centrality 
dependence of localized charged-neutral multiplicity 
fluctuations. It is based on an analysis of event-by- 
event fluctuation in the relative number of charged 
particles and photons detected within the common 
acceptance of the photon and charged particle mul- 
tiplicity detectors of the WA98 experiment . 

The paper is organized in the following manner: 
In the next section we describe the detectors used 
for the present analysis, the centrality selection cri- 
teria, the data reduction, and simulation. Section III 
deals with the analysis techniques, where two analy- 
sis methods are presented, one based on the correla- 
tion of photons and charged particles, and the other 
based on a discrete wavelet transformation analysis. 
In section IV, we present in detail the construction of 
mixed events used for this study. Section V discusses 
the ability of the mixed events to probe specific fluc- 
tuations. Final results and discussion are given in 
section VI. A summary is presented in section VII. 



II. EXPERIMENTAL SETUP AND DATA 
REDUCTION 

In the WA98 experiment at CERN [l^ , the main 
emphasis has been on high precision, simultaneous 
detection of both hadrons and photons. The experi- 
mental setup consisted of large acceptance hadron 
and photon spectrometers, detectors for charged 
particle and photon multiplicity measurements, and 
calorimeters for transverse and forward energy mea- 
surements. The present study makes use of the data 
from the photon multiplicity detector (PMD), the 
silicon pad multiplicity detector (SPMD) and the 
mid-rapidity calorimeter (MIRAC). 



A. Centrality Selection 

The centrality of the interaction was determined 
from the total transverse energy {Et) measured by 
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the mid-rapidity calorimeter (MIRAC) The 
MIRAC measures both the transverse electromag- 
netic and hadronic (iJ^"'^) energies in the in- 
terval 3.5 < 7^ < 5.5 with a resolution of 17.9%/a/E 
and 4:6.1%/ y/E, respectively, where E is expressed 
in GeV. The centrality, or impact parameter of the 
collision, has a strong correlation with the amount 
of i?T produced. Events with large Et production 
correspond to the most central, small impact param- 
eter, collisions pl| . 

The centralities are expressed as fractions of the 
minimum bias cross section as a function of the mea- 
sured total E't- For the present analysis we have 
used data selections in four centrality bins, the top 
5% (henceforth referred to as centrality-1), 5% - 
10% (centrality-2), 15% - 30% (centrality-3), and 
45% — 55% (centrality-4) of the minimum bias cross 
section. The minimum bias distribution of the total 
Et is shown in Fig. ^ The centrality bins used in 
this analysis are marked in the figure. 



B. Photon Multiplicity Detector 

The photon multiplicity is measured using the 
preshower photon multiplicity detector (PMD) lo- 
cated at a distance of 21.5 meters from the tar- 
get. The PMD consists of 3 radiation lengths (^o) 
thick lead converter plates in front of an array of 
square scintillator pads of four sizes, varying from 
15 mmx 15 mm to 25 mmx25 mm, placed in 28 box 
modules. Each box module consists of a matrix of 
38x50 scintillator pads read out using an image in- 
tensifier plus charged coupled device (CCD) cam- 
era system. The scintillation light is transmitted to 
the readout device via a short wavelength shifting 
fiber spliced to a long extra-mural absorber (EMA) 
coated clear fiber. The total light amplification of 
the readout system is ^40000. Digitization of the 
CCD pixel charge is done by a set of custom built 
fast-bus modules employing an 8 bit 20 MHz Flash 
ADC system. Details of the design and characteris- 
tics of the PMD may be found in Ref. (2^^J. The 
results presented here make use of the data from the 
central 22 box modules covering the pseudo-rapidity 
range of 2.9 < i] < 4.2. The clusters of hit pads, 
having total ADC content above a hadron rejection 
threshold are identified as photon-like, the multiplic- 
ity being denoted by Nj-nj^g. If the number of inci- 
dent photons is denoted by iV™'^ and the number of 
photons detected above the hadron rejection thresh- 
old as A^*'*, then the photon counting efficiency (e^) 
and purity of photon sample (/p) are defined as, 
= iV^V^r" and fp = Af^V^7-iikc respec- 
tively. These are estimated from detector simula- 
tions |2^ I23 . The photon counting efficiencies for 
the central to peripheral cases varies from 68% to 



73%. The purity of the photon sample in the two 
cases varies from 65% to 54% ^23, '23. The accep- 
tance in terms of transverse momentum (pr) extends 
down to about 30 M eV/c, however the PMD energy 
resolution is not sufficient for particle-by-particle pt 
measurement. 



C. Silicon Pad Multiplicity Detector 

The charged particle multiplicity (Ach) is mea- 
sured using the circular Silicon Pad Multiplicity De- 
tector (SPMD) located 32.8 cm from the target 
and having full azimuthal coverage in the region 
2.35 < 77 < 3.75, corresponding to the central ra- 
pidity region at SPS energies (where rjcMS — 2.92). 
The detector consists of four overlapping quadrants, 
each fabricated from a single 300 /im thick silicon 
wafer. The active area of each quadrant is divided 
into 1012 pads forming 46 azimuthal wedges and 22 
radial bins with a pad size increasing with radius 
to provide equal size pseudo-rapidity bins. The effi- 
ciency for detection of a charged particle in the ac- 
tive area has been determined in test beam measure- 
ments to be better than 99%. Conversely, the detec- 
tor is transparent to high energy photons, since only 
about 0.2% are expected to interact in the silicon. 
During the data recording, 95% of the pads worked 
properly and are used in this analysis. Details of the 
characteristics of the SPMD can be found in Ref. 
fisl 01 • The SPMD does not provide pr measure- 
ment, but provides the multiplicity measurement 
integrated over transverse momentum {j>t) with a 
threshold which extends down to about 20 MeV/ c. 



D. Data Reduction 

The data presented here were taken during De- 
cember 1996 at the CERN SPS with the 158-A GeV 
Pb ion beam on a Pb target of thickness 213 /Ltm. 
The WA98 Goliath magnet was switched off during 
these runs. Events with beam pile-up, downstream 
interactions, and pile-up in the CCD camera system 
were rejected in the off-line analysis 0, . The 
data have been analyzed for the region of common 
ri (2.9 < 7] < 3.75) and (j) coverage of the SPMD 
charged particle and PMD photon multiplicity de- 
tectors. The A^-,.-iikc and the A'ch distributions for 
the four centrality bins are shown in Figs. |21 and 
121 The number of events analyzed, the mean num- 
ber of photons and charged particles along with the 
root mean square deviations are shown in the figures. 
The PMD and SPMD detectors provide momentum 
integrated multiplicity measurements with very low 
thresholds. Since pions from DCC domains are ex- 
pected to have small pr values, below the pion mass, 
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the momentum integration will dilute the DCC sig- 
nal. On the other hand, the large coverage of the 
PMD and SPMD are important to overcome the lim- 
itations of small number fluctuations. 

The various sources of systematic errors associ- 
ated with the A^-y_iike and iVch distributions have 
been inv estig ated and described in detail previ- 
ously m Eg. These include: 

• (a) The uncertainty in the energy calibration 
and the associated uncertainty in the energy 
threshold for hadron rejection in the PMD 
leads to an error in the efficiency for A^-y-ukc 
clusters. The nominal hadron rejection thresh- 
old was set at three times the minimum ion- 
izing particle (MIP) peak. The value of the 
threshold was reduced by 10% 23] in order to 
estimate the systematic error. The associated 
error in 7V-y_iiko is 2.5%. 

• (b) The error due to the effect of clustering 
of pad signals in the PMD is a major source 
of error in A'^-uke- This error is determined 
from GEANT "2^ simulation by comparing 
the number of known photon tracks on the 
PMD with the total number of reconstructed 
photon-like clusters. It is found that the num- 
ber of clusters exceeds the number of tracks by 
3% in the case of peripheral events and by 7% 
for high multiplicity central events. 

• (c) The error due to the variation in pad- 
to-pad gains of the scintillators in PMD was 
found to be less than 1%. 

• (d) The uncertainty in the A^ch determina- 
tion with the SPMD has been estimated to be 
about 4% 

• (e) The error due to the finite resolution 
in the measurement of the total transverse 
energy(£'T) in MIRAC |2^ translates into an 
uncertainty in the centrality selection. The ef- 
fect of this systematic error has been deter- 
mined by performing the analysis with varying 
centrality cut within the MIRAC resolution. 

The contribution of each of these various systematic 
errors to the final results are discussed in the follow- 
ing sections. 



E. Simulated Events 

Simulated events were generated using the 
VENUS 4.12 event generator ^ with the de- 
fault parameter values. The output was pro- 
cessed through a detector simulation package in the 



GEANT 3.21 nil framework. This simulation in- 
cludes the full WA98 experimental setup and in- 
cludes experimental effects such as photon conver- 
sions, downstream interactions, hadron backgrounds 
in the PMD, etc. which might dilute or enhance the 
observed fiuctuations. The effect of Landau fluc- 
tuations in the energy loss of charged particles in 
silicon was included in the SPMD simulation 15]. 
For the PMD simulation, the GEANT results in 
terms of energy deposition in pads were converted 
to the pad ADC values using the MeV-ADC cali- 
bration relation. After this the ADC distribution 
is convoluted with a Gaussian function of proper 
width taken from the readout resolution curve. If 
the energy deposition is less than 3 MIP, a Landau 
distribution is used for convolution. The details of 
the PMD simulations taking into account the detec- 
tor and readout effects can be found in Ref. p^ . 
The centrality selection with the simulated data has 
been made in an identical manner to the data, de- 
termined from the simulated total transverse energy 
in MIRAC. The minimum bias total Et distribu- 
tion predicted by VENUS is shown by the dashed 
histogram in Fig.^ The agreement with the data is 
seen to be quite reasonable. A total of 60K VENUS 
events with simulated detector response were gener- 
ated for the present study. These simulated events 
(henceforth referred to simply as VENUS events un- 
less otherwise specified) were then processed with 
the same analysis codes as used for the analysis of 
the experimental data. 



III. ANALYSIS TECHNIQUES 

Two different analysis methods have been used 
in the present study. In the first analysis method, 
the magnitude of the iV-y_iike versus A^ch multiplic- 
ity fiuctuations is obtained in decreasing phase space 
regions. The second method employed the discrete 
wavelet transformation technique to investigate the 
relative magnitude of the 7V-),_iike versus TVch fluc- 
tuations in adjacent phase space regions. The re- 
sults from these methods of analysis applied to data, 
VENUS and various sets of mixed events (to be dis- 
cussed later) are compared to draw proper conclu- 
sions. 



A. A'^^ versus A^ch correlations 

In order to search for localized fluctuations in the 
photon and charged particle multiplicities, which 
may have non-statistical origin, the correlation be- 
tween Af-y-iiko and A^ch is investigated at various 
scales in 0. 
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The event-by-event correlation between A^-y-iikc 
and Nch has been studied in various 0-intervals by 
dividing the entire (p-apace into 2, 4, 8, and 16 bins. 
The method of analysis is similar to that described 
in Refs. ^3, 0|. Fig. 2| shows scatter plot of the 
correlation between A^-y_iike and A'ch for the top cen- 
trality bin. The correlation plots for each interval 
size, including the case of the full interval with no 
segmentation, are shown. The distributions for the 
other three centrality bins are qualitatively similar. 
A common correlation axis (Z) has been obtained 
for the full distribution by fitting the A'^-ukc and 
Nch correlation with a second order polynomial. The 
correlation axis with fit parameters is shown in the 
figure. The distance of separation (Dz) between a 
data point and the correlation axis has been calcu- 
lated with the convention that Dz is positive for 
points below the Z-axis. The distribution of Dz 
represents the relative fluctuations of A^^-iikc and 
Nch from the correlation axis for any chosen (j) bin 
size. In order to compare the fluctuations for dif- 
ferent (j) bins on a similar footing, a scaled vari- 
able, Sz — Dz/s{Dz), is used where s{Dz) rep- 
resents the rms deviation of the Dz distribution for 
VENUS events analyzed in the same manner. The 
Dz distributions of data, mixed events and the sim- 
ulated events for a given centrality and cj) bin size 
are all scaled by the same s{Dz) corresponding to 
the VENUS events for the respective centrality and 
azimuthal bin size. The presence of events with lo- 
calized fluctuations in iV^_iikc and Nch, at a given (j) 
bin, is expected to result in a broader distribution 
of Sz compared to those for normal events. Com- 
paring the rms deviations of the Sz distributions 
of data, mixed events (to be discussed later), and 
VENUS events may allow to infer the presence of 
non-statistical localized fluctuations. 



For the present DWT analysis the full azimuthal 
region is divided into smaller bins in cf), the number 
of bins at a given scale j being 2^ . The input to 
the analysis is a spectrum of the sample function at 
the smallest bin in 4> corresponding to the highest 
resolution scale, jmax- In the present case the sample 
function is chosen to be the photon fraction, given 
by: 

= N^-ui.M/iNy-meW + Nchm (2) 

A multi-resolution analysis has been carried out 
using the D—A wavelet basis on the above sam- 
ple function starting with jmax = 5. It may be 
mentioned that there are several families of wavelet 
bases distinguished by the number of coefficients 
and the level of iteration; we have used the fre- 
quently employed D—A wavelet basis 32]. The out- 
put of the DWT consists of a set of wavelet or fa- 
ther function coefficients (EEC) at each scale, from 
j — l,...,{jmax — !)■ The coefficients obtained at 
a given scale, j, are derived from the distribution 
of the sample function at one higher scale, j + 1- 
The EECs quantify the bin-to-bin fluctuations in the 
sample function at that higher scale relative to the 
average behaviour. The presence of localized non- 
statistical fluctuations will increase the rms devia- 
tion of the distribution of EECs and may result in 
non-Gaussian tails .29,. i3ll |. The DWT technique 
as used in this analysis has been demonstrated in 
our earlier publication '17*1. Once again, comparing 
the rms deviations of the EEC distributions of data, 
mixed events, and VENUS events may allow to infer 
the presence of localized fluctuations. The utility of 
the mixed events and the response of the analysis 
technique is demonstrated in Section V. 



B. Multi-resolution DWT analysis 

A multi-resolution analysis using discrete wavelet 
transformations (DWT) [23 has been shown to be 
quite powerful in the search for localized domains of 
DCC ,29, 30, SiJ. The beauty of the DWT technique 
lies in its power to analyze a spectrum at different 
resolutions with the ability to identify fluctuations 
present at any scale. This method has been uti- 
lized very successfully in many flelds including im- 
age processing, data compression, turbulence, hu- 
man vision, radar, and earthquake prediction p3 |. 
It should be noted that the DWT analysis provides 
different information than the moment analysis of 
the previous section. It analyzes the event-by-event 
distribution in phi space to characterize the bin-to- 
bin fluctuations relative to the average behaviour. 



IV. CONSTRUCTION OF MIXED EVENTS 

It is possible to search for non-statistical fluctu- 
ations in the experimental data in a model inde- 
pendent way by comparison of the data with mixed 
events generated from the data itself. Eurthermore, 
it is necessary to isolate the various contributions 
to the fluctuations and to understand all detector 
related effects in the data. This has been done 
by generating different types of mixed events which 
preserve the global multiplicity correlation between 
A^7-iikc and Nch- Eluctuations in the ratio of N^ to 
A^ch can arise due to ffuctuations in N-y only, fluc- 
tuations in A'ch only, or fluctuation in both N-y and 
A'ch- Furthermore, the ffuctuations in N^ and A'ch 
may be correlated event-by-event, as nominally ex- 
pected in the case of DCC formation. Each of these 
possibilities is investigated through the construction 
of four different kinds of mixed events. The method 
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of construction of these mixed events and the type 
of fluctuations they probe are described next. 



A. Maximally Mixed Events 

The first set of mixed events, referred to as Ml 
events, are constructed to remove afl correlations to 
the greatest extent possible to provide a baseline for 
comparison to the real event data. They were gener- 
ated by mixing hits in both the photon and charged 
particle detectors separately but still satisfying the 
global A'^-iikc-^ch correlation of the real event in 
the full acceptance. This means that on an event-by- 
event basis the total photon multiplicity and charged 
particle multiplicity of the mixed event were identi- 
cal to those of the real event to which it corresponds. 
Thus, the scatter plot of the mixed events is identi- 
cal to the real events for the single bin case shown in 
Fig. 0] The idea is to constrain the mixed events to 
be identical to real events globally and then compare 
them to the real data in localized regions of phase 
space to search for indications of non-statistical lo- 
calized fluctuations in the data. 

The Ml type of mixed events were constructed 
from the pool of all photon-like and charged particle 
hits, in which the hit position (rj and (f>) and event 
information (event number and total multiplicity) 
was kept for both the photon-like and charged par- 
ticle hits. For a given real event measured to have 
multiplicities Nj-m^e and Nch in the full acceptance 
region, a mixed event was constructed by randomly 
selecting iV^_iiko photon-like hits from the pool of 
photon-like hits and iVch charged particle hits from 
the pool of charged particle hits. This procedure 
was repeated for each real event. Care was taken 
such that no two hits from the same real event were 
used in the construction of a mixed event. Also, for 
the mixed events, hits within either the SPMD or 
PMD detector were not allowed to lie within the two- 
track resolution of that detector. In brief, the Ml 
mixed events randomly distribute the hits in each 
individual detector but keep the global correlation 
between the iV^-iikc and A^ch multiplicity. They pro- 
vide a maximally randomized sample of PMD and 
SPMD hits. Comparisons of such mixed events to 
real events will be most sensitive to the presence of 
localized fluctuations. However, in themselves they 
would not isolate the source of fluctuations as being 
due to and/or Nch, or correlations between 
and Nch- 



B. Minimally Mixed Events 

A second type of mixed events, referred to as 
M2 mixed events, were constructed to investigate 



the presence of correlated event-by-event fluctua- 
tions between Nj and Nch- These mixed events 
had a minimal amount of randomization since they 
were generated by mixing the photon hits taken un- 
altered from one event with the charged particle hits 
taken unmodifled from another event. As with the 
Ml mixed events, the global A^_iiko-^ch multiplicity 
correlation was maintained to be exactly the same as 
the data. To construct such mixed events, for a given 
real event measured to have multiplicities A^^uke 
and Nch in the full acceptance region, a mixed event 
was constructed by keeping the PMD A^_iike por- 
tion of the event intact, but combining it with the 
unaltered SPMD portion of a different randomly se- 
lected event, but constrained to have almost the 
same charged particle multiplicity Nch- This pro- 
cedure was repeated for each real event. In brief, 
this type of mixed event keeps the event-by-event 
hit structure in each detector identical to that in real 
events. Thus, such mixed events keep the individ- 
ual localized fluctuations present in N^ or Ach, but 
remove the event- by- event localized correlated fluc- 
tuations between them. Comparison of such mixed 
events to real events may reveal the presence of cor- 
related localized fluctuations between A^ and Ach. 



C. Partially Mixed Events 

Intermediate between the Ml and M2 types of 
mixed events are a third and fourth type of par- 
tially mixed event, referred to as M3-7 and M3-ch 
mixed events. These were constructed to provide in- 
formation regarding the contribution to the localized 
fluctuations in the A^ to Ach ratio from the individ- 
ual observables (A^ and Ach). They were generated 
from real events by mixing hits in one of the de- 
tectors (following the procedure for construction of 
Ml mixed events) and keeping the hit structure of 
the event in the other detector intact. M3-7 mixed 
events correspond to the case where the hits within 
the photon detector are unaltered while the hits in 
the charged particle detector are mixed. Similarly in 
M3-ch mixed events the hits in the charged particle 
detector were unaltered and the hits in the photon 
detector were mixed. In each type of mixed event the 
global A^_hke~-^ch correlation is maintained as in 
the real event. The two track resolution in the detec- 
tors where the hits are mixed is kept identical to that 
in real events. The total number of mixed events is 
the same as the number of real events. Comparison 
of such mixed events to real events and the other 
types of mixed events will reveal the presence of lo- 
calized fluctuations in A-y or Ach separately. 

A summary of the different sources of fluctuations 
in the ratio of A^ to Ach probed by each of the types 
of mixed events is given in Table 1 . 
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TABLE I: Type of fluctuations preserved by various 
mixed events 



Fluctuation 




Mixed 


Event 






Ml 


M2 


M3-ch 


M3-7 


Ai'-y-only 


No 


Yes 


No 


Yes 


Neil — only 


No 


Yes 


Yes 


No 


correlated N^-Nch 


No 


No 


No 


No 



V. DEMONSTRATION OF ANALYSIS 
METHOD 

In this section we wish to demonstrate the anal- 
ysis method and illustrate how the relationship of 
the measured result to that obtained with the var- 
ious mixed events can be used to provide an es- 
sentially model-independent signature of DCC for- 
mation. To demonstrate the analysis method and 
the potential to observe DCC event formation, we 
have applied the DWT analysis to a simple DCC- 
like model. The analysis is applied to "real" DCC 
events from the model as well as the various types of 
mixed events described in the previous section con- 
structed from the model DCC events. Since event 
generators which include DCC formation do not ex- 
ist, we have implemented a simple DCC model in 
which localized non-statistical N^-Nch fluctuations 
have been introduced by modification of the output 
of the VENUS event generator. To implement the 
fluctuations, the final state pions within a localized 
rj-fj) region from VENUS are interchanged pairwise 
(tt+tt" ^ 7r°7r") according to the DCC probability 
distribution, P{f) = 1/2^/f. The fiuctuations were 
generated over a localized region of 77 = 3 — 4 and a 
Aip interval of 90*^. The tt'^'s were then allowed to de- 
cay. The resulting events were then passed through 
the WA98 detector response simulation. The DCC 
events in the simple model used here give rise to 
an anti-correlation between and Nch- It also 
results in non-random fluctuations in both and 
Nch individually. Since the probability to produce 
events with localized charged-neutral fluctuations is 
unknown, ensembles of events, here referred to as 
"nDCC events" , were produced as a mixture of nor- 
mal VENUS events and events with localized fluctu- 
ations. The fraction of events with localized DCC- 
like fluctuations in each nDCC sample was varied as 
a parameter to be studied. 

By using VENUS events as the basis to introduce 
the DCC effect, it is insured that the general fea- 
tures of the event, including the multiplicity, com- 
position, and momenta of the produced particles, as 
well as correlations in the particle multiplicities due 
to impact parameter variation, are reasonably well 
described. Also the GEANT simulation of the de- 
tector response to these events insures that other 



effects such as photon conversions and the response 
of the PMD to hadrons, which might affect the ob- 
served multiplicities and the observed correlations, 
are also taken into account. While the assumption 
of a 100% DCC contribution over an interval of fixed 
size in 77 — (/) is certainly a gross simplification, it pro- 
vides a well-defined reference to gauge the potential 
for DCC observation. Other assumptions, such as 
a varying fraction, which might also be momentum 
dependent (since DCCs are expected to be a low pT 
phenomena) , and different or varying sizes might be 
more realistic. However, without clear theoretical 
guidance we have chosen this very simple model as 
a reasonable and well-defined reference. 

The DWT analysis was carried out on an ensem- 
ble of nDCC events and their corresponding mixed 
event sets created from each set of nDCC events. 
By ensemble of nDCC events we mean sets of events 
having different percentages of events with localized 
fiuctuations. The percentage varied from zero, that 
is, normal VENUS events with no localized fiuctua- 
tions, to an event set where all events had localized 
charged-neutral fiuctuations. For every nDCC event 
set the four sets of mixed events (M1,M2, M3-7, and 
M3-ch) were constructed from the nDCC events and 
analyzed. The rms deviations of the EEC distribu- 
tions for each set of nDCC events and corresponding 
mixed event sets were obtained. The results for scale 
j = I are plotted in Fig. as a function of the per- 
centage of localized DCC-type events in the nDCC 
event set. 

It is seen that the rms deviations of the FFC 
distributions of the nDCC events increase as the 
percentage of events with localized fiuctuations in 
N.y-Nch increases. This is the expected behaviour 
and demonstrates the linear response of the DWT 
to the DCC events when the frequency of events 
with fiuctuations increases. On the other hand, the 
rms deviations of the FFC distributions of the Ml 
type mixed events created from the nDCC events 
are found to be independent of the percentage of 
events having localized fluctuations. This is also the 
expected behaviour and demonstrates that the Ml 
mixed events can be used as a baseline from which 
to deduce the presence of fiuctuations in a model- 
independent manner. However, the deviation of the 
real events from the Ml mixed events does not in- 
form about the relative contributions of the individ- 
ual A^^ and N^h fiuctuations. The rms deviations 
of the FFC distributions of the M3 mixed events 
are found to be intermediate to those obtained for 
nDCC events and Ml mixed events. They indicate 
the separate contributions of the or Nch fiuctu- 
ations alone to the ratio. The rms deviations of the 
M2 mixed events are higher than those of Ml and 
M3 mixed events. That is because the M2 mixed 
events keep the separate contributions of both the 
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Nj and Nch fluctuations. However, the rms devia- 
tions of the M2 mixed events are consistently below 
those for the nDCC events. This is because the M2 
events randomize the correlations between Nj and 
Nch- The difference between the M2 mixed events 
and the nDCC events indicates the presence of the 
DCC-like correlated N^-Nch fluctuations. The rel- 
ative pattern of rms values for the real events and 
the various mixed events constructed from those real 
events seen in Fig. [S] provides a rather unambiguous 
model-independent signature for DCC-like fluctua- 
tions. Similar relative pattern of the rms deviations 
of the Sz distributions for mixed events and simu- 
lated events were also observed for N^-Nch correla- 
tion analysis. In particular, the observation of fluc- 
tuations in real events that are greater than the M2 
mixed events would constitute what might be called 
"smoking gun" evidence for DCC formation. Con- 
versely, the lack of a difference between real events 
and M2 mixed events would indicate the lack of 
DCC-like correlated charged-neutral ffuctuations. 



Figure |31 also demonstrates an effect which must 
be taken into account when comparing the measured 
result to the mixed events. For nDCC events with 
vanishing fraction of events with fluctuations, which 
is to say for normal VENUS events, it is seen that 
the rms deviations of the FFC distributions for all 
types of mixed events are higher than those for the 
nDCC events. In the VENUS simulations this is 
due to the presence of correlations between Nch and 
N^^-iike- These are primarily due to residual impact 
parameter correlations (see Fig. 0)) as well as due to 
the charged particle contamination in the N^-uke 
data sample whereby the charged particles register 
in both the PMD and SPMD (see Section II.C) _33j. 
These correlations are removed by the event mixing 
procedure which results in a larger rms deviations 
for the mixed events. The real data is presumed 
to have similar residual correlations as observed in 
the VENUS simulations. In order to correct for 
the effects of these non-DCC correlations, all mixed 
event rms values constructed from real events have 
been rescaled by the percentage difference between 
the rms deviations of the VENUS distributions and 
those of the corresponding VENUS mixed events, as 
also discussed in Ref. [T^ . 



For nDCC events with larger percentages of DCC- 
like events, the anti-correlation between Nch and 
Nj-uke overcomes the correlations between Nch and 
Nj-\ike and hence the rms deviations of the FFC 
distributions of nDCC events become greater than 
those of the mixed events, despite the other correla- 
tion effects. 



VI. RESULTS AND DISCUSSION 

In the analysis of experimental data, the results 
from the measured data are compared with simu- 
lated and mixed events. Below we discuss the results 
obtained from such a comparison using two different 
analysis methods discussed earlier. 



A. N-, versus A'^ch correlation results 

The Sz distributions calculated for different 4> bin 
sizes are shown in Fig.Elfor data, Ml, and VENUS 
events, for the four different centrality selections. 
The distributions for the other types of mixed events 
are not shown for clarity of presentation. The small 
differences in the Sz distributions have been quanti- 
fied in terms of the corresponding rms deviations, of 
these distributions shown in Fig. [7| The statistical 
errors on the values are small and are within the size 
of the symbols. The bars represent statistical and 
systematic errors added in quadrature. The various 
sources of systematic error have been discussed in 
an earlier section. An additional systematic error 
from the fit errors associated with the determina- 
tion of the correlation axis (Z) is also included. In 
general, the width of the Sz distribution increases 
from the most central to less central event selections 
and decreases with decreasing bin size. 

The rms deviations of the Sz distribution for the 
VENUS simulated events are 1 (by definition) for 
all centrality and all bins in azimuth by definition of 
Sz and are significantly different from the measured 
results. This is primarily because the global particle 
multiplicity as predicted by VENUS and those mea- 
sured in the experiment within the coverage of the 
detectors are not same. 

As seen in Fig. [7| the widths of the Sz distri- 
butions for mixed events closely follow those of the 
data. The mixed events have been constructed such 
that the global A'^-uko vs. Nch multiplicity correla- 
tions are maintained. Therefore the rms deviations 
of the data and the mixed event reference are the 
same in the first </> bin of Fig. [T] by construction. 
Some correlations between iVch and A'^-ukc are ex- 
pected, mostly as a result of the charged particle 
contamination in the iV^-ukc data sample, but are 
removed by the event mixing procedure and thereby 
result in a small difference between the real and 
mixed events, as seen in the analysis of the VENUS 
events, discussed earlier. All of the mixed event 
Sz distribution rms values ( Fig. 0) have therefore 
been rescaled by the percentage difference between 
the rms deviations of the VENUS distributions and 
those of the corresponding VENUS mixed events for 
each centrality class (Ref. [l^l. 
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The rms deviations of Sz distributions of the M2 
mixed events are found to agree with those of the ex- 
perimental data within errors for aU four centrahty 
classes and for all azimuthal bin sizes. This indi- 
cates the absence of event-by-event localized corre- 
lated fluctuations in A^^-uko and iVch, such as would 
be expected for DCC-like fluctuations. On the other 
hand, the rms deviations of the Ml mixed events are 
found to be systematically lower than those of the 
data for 2, 4, and 8 bins in for centrality bins 1, 
2, and 3. The results for both types of M3 mixed 
events are found to be intermediate between those 
of the data and the Ml mixed events. The results 
indicate the presence of localized fluctuations in the 
data in both the photon and charged particle mul- 
tiplicities. For the case of the most peripheral cen- 
trality selection (centrality-4) , the rms deviations of 
the Sz distributions of data and the various mixed 
events are found to be in close agreement to each 
other within the quoted errors. 



B. Multi-resolution DWT analysis results 

The FFC distributions, at scales j=l to 4, cor- 
responding to 4 to 32 bins in azimuthal angle, are 
shown in Fig. |H1 for data, Ml, and VENUS events 
for the four centrality classes. The results for other 
types of mixed events are not shown for clarity of 
presentation. The widths of the FFC distributions 
are found to increase in going from the most central 
to most peripheral centrality class. The rms devi- 
ations of these FFC distributions are summarized 
in Fig. El Similar to the case for the Sz distri- 
butions discussed above, the rms deviations of the 
mixed events have been rescaled by the percentage 
difference between the rms deviations of the VENUS 
FFC distributions and those of the VENUS mixed 
events for each centrality class. The statistical er- 
rors are small and are within the size of the symbols. 
The bars represent statistical and systematic errors 
added in quadrature. 

The rms deviations of the FFC distributions for 
the data, VENUS, and mixed events are found to be 
close to each other (within quoted errors) for the case 
of 32 bins in (f) for all of the four centrality classes. 
The rms deviations for the FFC distribution of M2 
mixed events are found to closely follow those of the 
data for all centrality classes and all bins in 0, while 
the rms deviations for the M3 mixed events lie be- 
tween those of the data and Ml mixed events. These 
results are consistent with those obtained from the 
analysis of the Sz distributions. These observations 
indicate the absence of event-by-event localized cor- 
related fluctuations (DCC-like) between N^^wi^a and 
A^ch- They also suggest the presence of localized 
fluctuations in both photons and charged particle 



multiplicities for intermediate bin sizes in azimuth. 
The rms values of the FFC distributions for VENUS 
events are close to those of the Ml mixed events 
for centrality classes 1, 2, and 3. However they 
are slightly higher for the most peripheral central- 
ity class (centrality-4). 



C. Discussion 

The results from the two independent methods 
of analysis are consistent and indicate the absence 
of event-by-event correlated DCC-like fluctuations 
in the photon and charged particle multiplicities. 
However they do suggest the presence of uncorre- 
lated fluctuations in both the photon and charged 
particle multiplicities for intermediate bin sizes in 
(f). The data has been compared to various kinds 
of mixed events and to simulated events which take 
into account many detector related effects. Still it 
is worthwhile to explore the extent of other possible 
experimental effects which might affect the observed 
rms deviations of the Sz and FFC distributions. As 
discussed extensively in Refs. [TEl l2^ l2^. care has 
been taken during the data taking and during the 
data processing to closely monitor the performance 
of the PMD and SPMD. The detector uniformity of 
the PMD was studied in detail by using the mini- 
mum ionizing particle (MIP) signal from the data 
for all pads of the 22 boxes. The fluctuations of the 
pad-to-pad relative gains were approximately 10%. 
The gain corrections were made for each pad. Cor- 
rections for gain variations during the data taking 
period were made periodically for both the PMD and 
SPMD. This reduces the possibility of abrupt gain 
or threshold changes during the run period. Events 
with obvious detector readout effects, such as miss- 
ing or dead regions, were carefully removed from the 
data sample. It should be recalled that most detec- 
tor effects are reflected in the mixed events. 

The effect of local fluctuations in the perfor- 
mance of the PMD has been studied using simulated 
VENUS events. In one test, the gain of a group of 
pads corresponding to one or more PMD cameras 
was randomly varied. A 30% change of gain in one 
camera for all events resulted in an increase of 0.7% 
in the rms deviations. Changing the gains of three 
cameras by 30% for all events resulted in an increase 
of the rms deviations by 1.7%. Because the cam- 
era gains were closely monitored on-line and during 
the processing of the data this is considered to be 
a highly unlikely scenario. Still, these changes are 
within the quoted errors of Figs. Eland l^which indi- 
cates that local gain and threshold variations would 
not account for the observed differences between the 
data and the Ml mixed events. 



10 



Several checks were performed to verify the qual- 
ity of the data obtained with the SPMD. One of the 
differences between the previous analysis |0 and 
the present one is that, while in the previous case 
the analysis was performed using the total charged 
particle multiplicity of the detector deduced from 
the magnitude of the measured SPMD signals, here 
we simply used the total number of hit SPMD pads. 
The reason for using hit pads is that the correc- 
tion in going from deposited charge to hits for each 
event and small r]-(j} segments is non-trivial. Also, 
the effects of two-track resolution and possible shifts 
of the beam position on the target during the spill 
were studied in detail. However all of these pro- 
duced small effects which could not account for the 
differences observed between the data and Ml mixed 
events. 



D. Strength of localized fluctuations 

In order to quantify the strength of the iV^_iiko 
and A'ch fluctuations for various bins in (p and for 
different centrality classes, we define a quantity x 
as: 

Sl 

where si and s correspond to the rms deviations 
of the FFC distributions of the Ml mixed events 
and real data, respectively. The results are shown 
in Fig. ^] as a function of the number of bins in (j> 
for the four different centrality classes. Qualitatively 
similar results are obtained when x is calculated us- 
ing the rms deviations of the Sz distributions. The 
shaded portion indicates the region of x where s is 
one a greater than the rms deviation FFC distri- 
butions for Ml events, where a is the total error 
on the Ml event rms deviation. It represents the 
limit above which a signal is detectable. Since x is 
calculated from the rms deviations of the FFC dis- 
tributions for data and Ml mixed events it gives the 
combined strength of localized fluctuations in both 
the photon and charged particle multiplicities. We 
do not present x values calculated using M3-type 
mixed events. However, it is clear from the rms de- 
viation figures (Fig.[7|andO that both photons and 
charged particles contribute to the observed fluctu- 
ations. The result shows that the strength of the 
fluctuations decreases as the number of bins in (j) 
increases, with a strength which decreases to below 
detectable level (within the quoted errors) for 16 and 
32 bins. There is an indication that the strength of 
the signal decreases with decreasing centrality for 4 
and 8 bins in azimuthal angle, although the tendency 
is not very strong. 



E. Upper limit on DCC production 

It has been shown that the rms deviations of the 
Sz and FFC distributions for data are very close 
to those of the M2 mixed events, within the quoted 
errors. If the DCC-like correlated fluctuations in 
iV^_iikc vs. Nch were large, the rms deviations (Fig.[7| 
and Fig.O of data would have been larger compared 
to those of the M2 mixed events. Since this is not the 
case, we may extract an upper limit on the produc- 
tion of DCCs at the 90% confidence level following 
the standard procedure as discussed in Ref. 34] . 

The errors are assumed to have Gaussian distri- 
bution, although they are asymmetric. The larger of 
the asymmetric errors is conservatively used for the 
limit calculation. The 90% C.L upper limit contour 
has been calculated as x + 1.28e;^, where x is calcu- 
lated using Eq.®. Here si and s correspond to the 
rms deviations of the FFC (or Sz) distributions for 
M2 mixed events and real data, respectively, and 
is the error in x from the FFC (or Sz) analysis. If 
the difference between the rms deviation of the FFC 
distributions for M2 mixed events and real data is 
negative, we take the value of x to be zero. It may 
be mentioned that for the calculation of the upper 
limits we have assumed that the total difference in 
the rms values of data and M2 mixed events is due 
to DCC-like fluctuations only. 

To relate the measured upper limit on the size 
of the fluctuations to a limit on DCC domain size 
and frequency of occurrence we proceed as follows: 
Within the context of the simple simulated DCC 
model described earlier, we obtain the rms devia- 
tions of the FFC distributions with various domain 
sizes in azimuthal angle (15", 30", • • • 180") and 
for each domain size also for different frequency of 
occurrence of DCC (0% to 100%). The M2 mixed 
events are then constructed for each of these sets of 
simulated events. For each set of DCC-type events 
of a given domain size and frequency of occurrence, 
the value of x is calculated using Eq. 13 from the dif- 
ference in rms deviations of the FFC distribution of 
the DCC event (s) and its corresponding M2 mixed 
event distribution (sl). The upper limit is set at that 
value of frequency of occurrence for a fixed DCC do- 
main size at which the x value from the DCC model 
matches with that of the x + 1.286^^ upper limit from 
the experimental data. This is used to set the upper 
limit contour in terms of domain size and frequency 
of occurrence of the DCC. The results for centrality 
classes 1 and 2 are shown in Fig. It may be men- 
tioned that the upper limit contour set by a similar 
analysis of the rms values of the Sz distributions is 
very similar to that from the FFC analysis. Also, 
from Fig. ^| it is seen that the total fluctuation is 
similar or weaker for centrality classes 3 and 4, hence 
the upper limit for these two classes would be similar 
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or weaker than those shown m Fig. II II for centrahty 
classes 1 and 2. 



VII. SUMMARY 

A detailed event-by-event analysis of the rj — 
4> phase space distributions of the multiplicity of 
charged particles and photons in Pb+Pb collisions 
at 158-^ GeV has been carried out using two dif- 
ferent analysis methods for four different centrality 
classes. The results from the two analysis methods 
were found to be consistent with each other. The 
first analysis method studied the magnitude of the 
-^7 -like versus A^ch multiplicity fluctuations in de- 
creasing phase space regions. The second analysis 
employed the discrete wavelet transformation tech- 
nique to investigate the relative magnitude of the 
A^7-iikc versus iVch fluctuations in adjacent phase 
space regions. The results were compared to pure 
VENUS-I-GEANT simulation events and to various 
types of mixed events to search for and identify the 
source of non-statistical fluctuations. Both analy- 
sis methods indicated fluctuations beyond those ob- 
served in simulated and fully mixed events for (j) 
intervals of greater than 45° and which increased 
weakly in strength with increasing centrality. The 
additional fluctuations were found to be due to un- 
correlated fluctuations in both Af-^-ukc and A^ch- No 
significant correlated fluctuations in A'^-ukc versus 
A^ch, a likely signature of formation of disoriented 
chiral condensates, were observed in all of the four 
centrality classes studied. Using the results from the 
data, mixed events, and within the limitations of a 
simple model of DCC formation, an upper limit on 
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FIG. 1; The total Et distribution (solid dots) measured in 3.5 < j? < 5.5 for Pb+Pb collision at 158- A GeV/c. The 
total Et distribution as obtained from VENUS is also shown as dashed histogram. The Et values corresponding to 
the different centrality bins used in the analysis are shown, (a) centrality-1 (0 — 5%), (b) centrality-2 (5 — 10%), (c) 
centrality-3 (15 — 30%), (d) ccntrality-4 (45 — 55%) of the minimum bias cross section as determined by selection on 
the measured transverse energy distribution. 




FIG. 2: The A'^.-iike multiplicity distributions for the four centrality selections for Pb+Pb collision at 158- A GeV/c. 
(a) centrality-1 (0 - 5%), (b) centrality-2 (5 - 10%), (c) centrality-3 (15 - 30%), (d) centrality-4 (45 - 55%). 





FIG. 3: The A'^ch multiplicity distributions for the four centrality selections for Pb+Pb collision at 158- A GeV/c. 
centrality-1 (0 - 5%), (b) centrality-2 (5 - 10%), (c) centrality-3 (15 - 30%), (d) centrality-4 (45 - 55%). 
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FIG. 4: The event-by-cvcnt correlation between A'ch and A'^^-iikc for centrality bin=l Overlaid on the plot is the 
common correlation axis (Z-axis) obtained for the full distribution by fitting the A''^_iikc and A'^ch correlation with a 
second order polynomial. The values in the upper left are the coefficients of the polynomial of the form A'ch = a + 
biV-^-iike+cA/^-iike- For the other three centrality classes, the plots are qualitatively similar. 
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FIG. 5: The rms deviations of the FFC distributions for 
simulated VENUS events containing a variable fraction 
of localized DCC- like events with DCC extent Ac/)_dcc = 
90°, as a function of that fraction. Results are also shown 
for various mixed events constructed from those events. 
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FIG. 6: The Sz distributions for data (solid circles), mixed events (solid histogram), and simulated events (dashed 
histogram) for the four centrality bins (a) centrality-1 (0 — 5%), (b) centrality-2 (5 — 10%), (c) centrality-3 (15 — 30%), 
(d) centrality-4 (45 - 55%). 
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FIG. 7: The RMS deviations of the Sz distributions for data, various mixed events, and simulated events for the 
four centrality bins (a) centrality-1 (0 — 5%), (b) centrality-2 (5 — 10%), (c) centrality-3 (15 — 30%), (d) centrality-4 
(45 - 55%). 
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FIG. 8: The FFC distributions for data (solid circles), mixed events (solid histogram), and simulation (dashed 
histogram) for the four centrality bins (a) centrality-1 (0 — 5%), (b) centrality-2 (5 — 10%), (c) centrality-3 (15 — 30%), 
(d) centrality-4 (45 - 55%). 
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FIG. 9: The RMS deviations of the FFC distributions for data, various mixed events, and simulation for the four 
centrality bins (a) centrality-1 (0 — 5%), (b) centrality-2 (5 — 10%), (c) centrality-3 (15 — 30%), (d) centrality-4 
(45 - 55%). 
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FIG. 10: The fluctuation strength parameter for the four ccntrahty classes. Ccntrahty-1 corresponds to the 5% most 
central, centrality-2 corresponds to 5 — 10%, centrality-3 corresponds to 15 — 30% and centrality-4 corresponds to 
45 — 55% of the minimum bias cross section as determined by selection on the measured transverse energy distribution. 
The error bars axe shown only on the centrality-1 selection for clarity of presentation. The errors are similar for the 
other centralities. The shaded portion represents the limit above which a signal is detectable (see text for details). 
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FIG. 11: The 90% confidence level upper limit on DCC production for central Pb+Pb collision at 158- A GcV/c, as 
a function of the DCC domain size in azimuthal angle within the context of a simple DCC model and the measured 
photon and charged particle nmltiplicities in the interval 2.9 < r] < 3.75. The solid line corresponds to data from the 
top 5% and dashed line to top 5 — 10% of the minimum bias cross section as determined by selection on the measured 
transverse energy distribution. 



